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Frame Dependence of Parton Cascade results 
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Frame dependence of parton cascade results is studied for different schemes of doing cascade sim- 
ulations. We show that different schemes do not always agree and results may have strong frame 
dependence. When the interaction range is on the order of mean free path and the collisions are 
done in the two parton center of mass frame, the results are not sensitive to the global frame 
that the collisions are ordered. Effects of different differential cross sections are also discussed. 



I. INTRODUCTION 

Cascade models have been widely used to simulate 
Rclativistic and Ultra-Relativistic Nucleus-Nucleus col- 
lisions JO. They can be used to simulate systems that 
are not in thermal equilibrium and they incooperate the 
effect of finite mean free path automatically. 

Causality violation is inherent in cascade simulations 
due to the geometrical interpretation of the cross sec- 
tion, i.e., particles scatter when their closest distance is 
smaller than (cr is the scattering cross section). One 
problem is that the information travels across at a 
time which is either the fixed time step or the mean free 
path. This may lead to a speed of information faster 
than the speed of light and may cause superluminous 
shock waves. Another problem is that different choices 
of doing collisions, i.e., different collision schemes, may 
lead to different collision orderings and hence different 
physical predictions. 

Many groups have studied the causality problem for 
energies below RHIC energy g] and attempts have been 
made to reduce the problem. For Ultra-relativistic heavy 
ion collisions at RHIC energies and beyond, causality vi- 
olation may be very serious due to the relatively small 
mean free path which may be on the order or even much 
less than the interaction range. The superluminous sig- 
nals have been studied Q in cascade simulations of pA 
collisions. In this paper, we study the effect of colli- 
sion ordering on macroscopic variables by comparing the 
results of different collision schemes. A parton cascade 
code, GPC (Generic Parton Cascade), we developed re- 
cently Q has been used in this study. 

First, we describe the initial conditions and different 
collision schemes. Then we compare the results from dif- 
ferent schemes with a discussion of the effects of different 



differential cross sections followed by the conclusions. 



II. INITIAL CONDITIONS AND COLLISION 
SCHEMES 

To study frame dependence of cascade results at RHIC 
energies, we prepare a system of gluons similar to the 
minijet gluon system that is going to be produced at 
RHIC. 4000 gluons are uniformly distributed in the —5 to 
5 space time rapidity range. Initially, they get local ther- 
mal equilibrium at temperature 500 MeV. They occupy 
a transverse disk of radius 5 fm at t — and are formed 
at longitudinal proper time To = 0.1 fm (i.e., the forma- 
tion time for a particular particle i is ti = tq x cosh(rji), 
Zi = t x sinh(r]i)). 

We will compare results in different frames of ordering 
the collisions for the following schemes of doing cascade: 

(a) Two gluons collide when their closest approach dis- 
tance in the two particle center of mass frame is smaller 
than (a is the scattering cross section) . The collision 
space point is chosen to be the midpoint of the two parti- 
cles in the two-body center of mass frame at their closest 
distance, the collisions are ordered in a global frame, ei- 
ther the collider lab frame or the target frame (they differ 
by 6 unit of rapidity boost); 

(b) differs from (a) in that the collision space point 
for a particle (not for a collision) is the position of the 
particle in the two-body center of mass frame. So, in 
the global frame, for a particular parton parton collision, 
each parton will have its own time of scattering. The 
collisions are ordered according to the average time of 
the two scattering times in the global frame; 

(c) differs from (b) in that the collisions are ordered 
according to the earlier time of the two scattering times; 
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(d) two partons collide when their closest distance is 
smaller than »/^ in the global frame. The collisions 
are ordered according to the collision time in the global 
frame. 

Different schemes specify different parton collision 
frame and different collision ordering time. 

The system propagates from one collision to the next 
in the global frame instead of from one time step to the 
next. 

Partons are all on mass shell. Rutherford cross section 
regulated by a screening mass, ^ = 2 ^°f p > is taken 
to generate the scattering angular distribution. In the 
above formula, t = (pi — P2) 2 in which p\ and P2 are 4 
momenta of one particle before and after the collision. 
The total cross section for parton collision is taken as: 
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We set = 

Numerically we divide the space into cells. We 
store only the next collision time, partner and next par- 
ticle in the same cell in an interaction list to save memory 
and memory manipulation. The interaction list is only 
updated locally, i.e., only the cell where the collision hap- 
pens and its neighboring cells arc updated. This reduces 
the number of checks and significantly increases the speed 
of the parton cascade code. We can get one 6,000 colli- 
sions event in around 10 minutes. Our results have been 
checked by taking away the cells to make sure they are 
consistent with the ordinary calculations without space 
divisions. At the present stage, there are only 2 to 2 
scatterings. We are working on including the radiation 
processes into the GPC. 



III. RESULTS FOR DIFFERENT SCHEMES 



In our case, R = 5/m, ^ = 400 and t ~ 0.2/m taken 
into account the spread of rapidity around the space time 
rapidity. We get n ~ 25//m 3 . 
The interaction range ~ — . 

V x ^ 

To have I > we need to have /i > 4.3 fm^ 1 m 

0.85 GeV . This is a condition for isolated 2 body col- 
lisions that is generally violated. Ideally we also need 
many particles within the Debye screening scale: 
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For/i = 3/m- 1 , N D « ^ w 4. 

We showed below even though the isolated 2 body col- 
lision condition is barely satisfied, in the situation with 
strong inside-outside correlations like the system evolving 
from the initial conditions we specified above, a sensible 
scheme can still produce reasonable results. 

Fig. 1, Fig. 2 and Fig. 3 give the ^, and 

distributions with cross section a = ^ fm 2 (a = 0.47 
and \x = 3 fm -1 ) respectively. In plotting the target 
frame ordering case, we make a shift so that the central 
rapidity always has rapidity 0. 

After evolution, ^P- is almost the same as the initial. 

^j^- has a drop from initial value ~ 385 GeV to around 
335 GeV . This shows the longitudinal work Q done by 
the effective pressure. The absolute value of slope of 
log curve is increasing which clearly shows the effective 
cooling of the system. The results do not depend on the 
frame we use to order the collisions. 

Scheme (b) and (c) give similar results. 

A close look at the average number of collisions per 
event shows that for scheme (a) we get ~ 6900 collisions, 
for (b) 6800 and for (c) 6200 when collision order is in 
the collider lab frame. In the target frame case, each gets 
several hundred less. 
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distribution of partons are 
closely related to the hadronic observables. In this 
section, we are going to study parton cascade induced 
changes of these distributions. First, we look at the case 
when the mean free path is on the order of the interac- 
tion range y^f; then the case when the mean free path 
is much smaller. Further, we study the dependence of 
the results on the differential cross section, i.e., the an- 
gular distribution of outgoing particles from a collision. 
Finally, a discussion on the number of total collisions vs. 
the number of non-causal collisions and effects of a small 
particle mass. 

The mean free path I ~ — in which n is the par- 
ton number density and the cross section a 
density n can be estimated through: 
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FIG. 1. 20 event averaged 4£L distribution for scheme (a) 
with a = 0.47 and /i = 3 fm~ . The dot-dashed curve is 
the initial distribution for reference; the solid curve is for the 
case that collisions are ordered in the collider center of mass 
frame; for the dashed curve, the collisions are ordered in the 
target frame. 
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FIG. 2. 20 event averaged distribution for scheme (a) 

with a — 0.47 and /i = 3 fmT 1 . The dot-dashed, solid and 
dashed curves are for the initial, collider frame ordering and 
target frame ordering respectively. 



in scheme (d), we get much more collisions even with a 
smaller cross section. In the collider frame ordering case, 
there are around 22, 000 collisions per event; in the target 
frame ordering case, we get 97, 000! 

The frame dependence of scheme (d) is due to the ra- 
pidity correlation introduced by taking the closest dis- 
tance in the global frame. It neglects the flow velocity 
and increases the density of the particles. This is also 
seen in the huge increase of number of collisions with 
smaller cross section comparing to the other schemes. 
Scheme (a), (b) and (c) take into account of the flow by 
looking at the closest distance in the 2 particle center 
of mass frame. This corresponds to the time in the 2 
particle center of mass frame when the 2 discs meet. 
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FIG. 3. 20 event averaged distribution for scheme (a) 
with a — 0.47 and /i = 3 fm~ . The dot-dashed, solid and 
dashed curves are for the initial, collider frame ordering and 
target frame ordering respectively. 

Fig. 4, Fig. 5 and Fig. 6 give the results for scheme (d). 
Strong frame dependence is observed in this case. From 
Fig. 4 and 5, we see that a boost of the system by 6 
unit of rapidity (the original rapidity goes to —6 before 
shifting back to 0) makes a peak which is shifted from 
the central along the boost by 2 unit of rapidity. We use 
a much smaller a = ^ fm 2 « lmb cross section because 



FIG. 4. 5 event averaged ^ distribution for scheme (d) 
with a = 0.47 and /i = 5 fmT 1 . The dot-dashed, solid and 
dashed curves are for the initial, collider frame ordering and 
target frame ordering respectively. 




FIG. 5. 5 event averaged distribution for scheme (d) 
with a = 0.47 and /i = 5 fmT 1 . The dot-dashed, solid and 
dashed curves are for the initial, collider frame ordering and 
target frame ordering respectively. 
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FIG. 7. 20 event averaged ^ distribution for scheme (a). 
The dot-dashed curve is for the initial distribution; the thick 
solid is for a = 0.47 and /i = 1 fm' 1 collider frame ordering; 
the thick dashed curve is for a = 0.47 and fi = 1 fm~ 
target frame ordering; the long dashed curve is for a = 1.41 
and jj, = 3 fm' 1 collider frame ordering; the dotted curve 
is for a — 1.41 and fi — 3 fm' 1 target frame ordering. The 
dot-dashed, thick dashed and dotted curves overlap; the thick 
solid and long dashed curves overlap. 



FIG. 6. 5 event averaged distribution for scheme (d) 

with a = 0.47 and /i = 5 fm' 1 . The dot-dashed, solid and 
dashed curves are for the initial, collider frame ordering and 
target frame ordering respectively. 

Fig. 7, Fig. 8 and Fig. 9 give results for a = ir fm 2 . 
Now the interaction range is 1/m and the mean free path 
is initially ^ ~ 0.01/m. For a — 0.47 and /i = 1 fm -1 , 
the dN/dy distribution has a dip in the middle for the 
collider frame ordering case, while the target frame order- 
ing case doesn't have one. The dEx/dy central rapidity 
plateau height is almost the same as the for a = 0.47 
and \i = 3 fm' 1 case. This is because of the effective 
screening mass is used to regulate the forward scattering 
cross section and larger cross section or more forward 
collisions does not lead to more collective work. This is 
clearly shown when we vary the differential cross section 
by change the a from 0.47 to 1.41 and /j, from 1 fm' 1 
to 3 fm' 1 . As shown in Fig. 8, the collider frame or- 
dering case has a clear drop at rapidity from —4 to —2 
and from 2 to 4. In the middle, there is a peak. The 
target frame ordering has a lower plateau. The dEx/dy 
is frame dependent now. The dN/dpx for a — 1.41 and 
fi = 3 fm' 1 has a lower temperature parameter than 
that of a = 0.47 and /i = 1 fm' 1 . This is because more 
work is been done when there are more relatively large 
angle scatterings and we get more cooling. Scheme (b) 
and (c) can not get frame independent results either. 
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FIG. 8. 20 event averaged ^jff- distribution for scheme (a). 
The dot-dashed curve is for the initial distribution; the thick 
solid is for a — 0.47 and fi — 1 fm collider frame ordering; 
the thick dashed curve is for a = 0.47 and /i = 1 fm' 1 target 
frame ordering; the long dashed curve is for a = 1.41 and 
/i = 3 fm' 1 collider frame ordering; the dotted curve is for 
a = 1.41 and fi = 3 fm' 1 target frame ordering. 
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FIG. 9. 20 event averaged -j^ distribution for scheme (a). 
The dot-dashed curve is for the initial distribution; the thick 
solid is for a = 0.47 and fi = 1 fm' 1 collider frame ordering; 
the thick dashed curve is for a = 0.47 and /i = 1 fm' 1 target 
frame ordering; the long dashed curve is for a = 1.41 and 
fi = 3 fm' 1 collider frame ordering; the dotted curve is for 
a = 1.41 and /i = 3 fm' 1 target frame ordering. 
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To study the effect of different differential cross sec- 
tions, we compare the results for a = 0.47 and fi = 
3 fm" 1 and a = 1.41 and ^ = 9 fm . Shown in 
Fig. 10 and Fig. 11 are dEr/dy and dN/dpT distribu- 
tions. dN / dy distribution is the same for these two cases 
and they are the same as the initial. From Fig. 10, a 
clear drop in the height of central rapidity plateau can 
be seen with more large angle scatterings. This shows 
that more work has been done and is also the reason 
that we get more cooling in dN/dpx distribution. Wc 
get 7, 000 collisions in the collider frame and 6, 700 in the 
target frame. Scheme (b) gives the similar results and 
they are all frame independent. 
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Another interesting thing to look at is the number of 
non-causal collisions vs. the number of total collisions. A 
non-causal collision here is defined to be a collision that 
the global ordering is not consistent with the ordering for 
a particular particle participating in the collision. This 
is better seen in scheme (b). If particle a's next collision 
according to global ordering is with i, but the collision 
time for particle a with i (i.e., the time a is going to 
change its momentum) is later than the collision time for 
a with j (even though according to global ordering, a j 
collision happens later), then we call this collision a non- 
causal collision. For scheme (b) in the collider lab frame, 
with a = 0.47 and fi = 3 fm" 1 , we get a total of 6,800 
collisions per event with 900 of them non-causal. With 
a larger cross section (a = 0.47 and \i = 1 fm" 1 ), as we 
expect, the non-causal to total ratio goes up. There are 
total 21,000 collisions with 10,400 of them non-causal. 
Boost to the target frame, the total number of collisions 
for a = 0.47 and /i = 3 fm" 1 is 6, 500 while the number 
of non-causal is 1, 000. The effects of putting in a mass 
for the particles are also studied. There are no significant 
changes in , ^ and and the number of total and 
non-causal collisions for massless, m — 0.01 GeV and 
to = 0.1 GeV. Fig. 12 shows the collisions rates. We see 
that the results are not sensitive to small masses of the 
particles. 



FIG. 10. 20 event averaged distribution for scheme 
(a). The dot-dashed curve is for the initial distribution; the 
dotted curve is for a — 0.47 and ^ = 3 fm" 1 the thick solid 
is for a = 1.41 and /i = 9 fm" 1 collider frame ordering; the 
thick dashed curve is for a = 1.41 and /i = 9 fm -1 target 
frame ordering. 



10000 




Pt (GeV) 

FIG. 11. 20 event averaged 4^- distribution for scheme 
(a). The dot-dashed curve is for the initial distribution; the 
dotted curve is for a — 0.47 and jj, = 3 fm" 1 the thick solid 
is for a = 1.41 and /i = 9 fm~ collider frame ordering; the 
thick dashed curve is for a = 1.41 and /i = 9 fm" 1 target 
frame ordering. 
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FIG. 12. 20 event averaged collision rate as a function of 
time for scheme (b). The upper 3 solid curves are for massless, 
m — O.OlGeV and m = O.lGeV total collision rate and the 
lower 3 dashed curves are for non-causal collision rate of these 
3 masses. 

What we describe here is the parallel ensemble method 
in which average is taken over n independent events. An 
alternative, the full ensemble method ||], is to increase 
the particle number by k (by dividing the particles into 
small pieces) and at the same time decrease the cross sec- 
tion by a factor of k. In the full ensemble method, the 
mean free path is the same as for the parallel ensemble 
method, but the cross section is smaller and hence iso- 
lated 2 body scattering condition can be satisfied. This 
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makes it possible to eliminate the causality violation in 
the limit of large k. But the full ensemble method will 
increase the computation time significantly due to the in- 
creasing number of checks required for the higher density 
system. We see from the above study that for reactions 
with strong inside outside correlations, the full ensemble 
method may not have to be used in practice to get reason- 
able results as long as a sensible scattering prescription 
is used. 



IV. CONCLUSIONS 
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After studying the frame dependence of parton cascade 
results, we see that different schemes may give quite dif- 
ferent results. Especially, the global frame collision and 
ordering scheme results have a strong frame dependence. 
Other 2 particle center of mass frame collision and global 
frame ordering schemes can give almost frame indepen- 
dent results even when the interaction range is on the 
order of the mean free path. As we expect, they fail 
when the interaction is much large than the mean free 
path. Different differential cross sections lead to different 
longitudinal work and give different dEx / dy and dN/ dpx 
results. Since radiation will increase the cross section and 
parton density, the quantum statistics will also change ef- 
fective cross sections, further study is necessary when we 
apply this to more realistic situations. 
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